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If general relativity (GR) describes the expansion of the Universe, the observed cosmic acceleration 
implies the existence of a 'dark energy'. However, while the Universe is on average homogeneous 
on large scales, it is inhomogeneous on smaller scales. While GR governs the dynamics of the 
inhomogeneous Universe, the averaged homogeneous Universe obeys modified Einstein equations. 
Can such modifications alone explain the acceleration? For a simple generic model with realistic ini- 
tial conditions, we show the answer to be 'no'. Averaging effects negligibly influence the cosmological 
dynamics. 



Cosmological observations have established that we 
live in an accelerating Universe Assuming the dy- 
namics of the expanding Universe to be described by gen- 
eral relativity (GR), the observed acceleration compels us 
to conclude the existence of a 'substance' with negative 
pressure named 'dark energy' 0. The standard model 
of particle physics offers no candidate for dark energy, 
nor does it provide any convincing clues as to its origin. 
While this situation might suggest new physics, it is still 
worthwhile to ask if a conventional explanation for the 
acceleration might yet be possible. 

It is true that the standard hot Big Bang model for 
a homogeneous and isotropic Universe described by the 
Friedmann-Lemaitre-Roberton- Walker (FLRW) solution 
of Einstein's equations agrees well with observations. 
However, the actual Universe is highly inhomogeneous 
and homogeneity is recovered only in an average sense 
on sufficiently large scales Q . Strictly speaking, only the 
inhomogeneous Universe obeys Einstein's equations, and 
recovering the homogeneous limit by averaging leads to 
corrections in the equations Could these correction 
terms (backreaction) be significant enough to account for 
dark energy? This question has attracted considerable 
recent attention 0, [a 0| • 

One might simply take the view that since the per- 
turbed FLRW model works well, backreaction from aver- 
aging these small perturbations cannot be large. There 
is however a subtle catch here. The evolution of inho- 
mogeneities depends on the background FLRW solution, 
which in turn depends on the backreaction. One cannot 
a priori rule out a runaway situation in which the back- 
reaction is fed by and reinforces a strong evolution of the 
inhomogeneities, perhaps leading to a dark energy like 
component at late times. 

A covariant averaging formalism for GR known as 
macroscopic gravity (MG) has been developed by Za- 
laletdinov [8(. In MG the connection on the inhomoge- 
neous manifold M. is averaged to obtain a connection 
on a new "averaged manifold" M in a rigorous manner. 
As a consequence, the averaging of Einstein's equations 
for M. leads to a set of corrected equations for M. with 
the corrections being quadratic in the underlying connec- 



tion. Of relevance for cosmology is the special case where 
averaging is performed on specific spatial slices; the back- 
reaction for this case was obtained in [9(. Further, if it 
can be self-consistently established that the perturbed 
FLRW model works well throughout the evolution of the 
Universe, then the backreaction must be evaluated for 
this specific perturbative case. This has been done in 
[icj |. by assuming the metric to have the form 
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with r\ the conformal time and jab the (flat) 3-space met- 
ric. In the present paper we work with the cosmic time 
r related to the conformal time by dr = a(r))drj, with 
appropriate conversions applied. The modified FLRW 
equations are given by 
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where p and p are the background energy density and 
pressure respectively and the combinations (V^ + S^) 
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and {V^ + + S^) are covariant scalars given by 



2({d T 4'f) + {{drV - dr^f) 

-((\7 A \7 B d T f3) (V A \7 B d T P)) 

&{d A 1}>d A 1p) + (d A (<p - 4;)d A & - vo> 

- ((V a VbV c /3)(V a V s V c /3))1 , 
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Here (A, B = 1, 2, 3) is the 3-space covariant deriva- 
tive and (3 solves V 2 /3 = <p - Zip with V 2 = ^ AB V A V B , 
with the condition that when ip = = if), (3 = 0. The 
angular brackets in © denote a spatial averaging defined 



(f)(r,x) = ^[ d 3 yf(r,y). 

V L Jv{x) 



(4) 



for any function f(r,x) where V(x) is the 3-dimensional 
averaging domain of "Eulerian" length scale L (comoving 
with the background) and volume Vl- This averaging 
operation is derived by considering a spatial averaging 
limit of the full MG averaging technique and working in 
a specific "volume preserving" gauge (in which the metric 
determinant depends only on cosmic time) with the form 



ds 2 = -(1 + 64>)dr 2 + a 2 (r)(l - 2^)da 



(5) 



where ip = tp/3 (see [lfl]). This procedure is thus dif- 
ferent from the spatial averaging of Buchert @, Q . Hav- 
ing defined the averaging in this gauge, one rewrites all 
quantities in terms of the conformal Newtonian gauge 
potentials ip and ip, and ([3]) show the final result. 

In this paper we work with the background metric in 
spherical coordinates, 



ds 2 bg - -dr 2 + a 2 (r) [dr 2 + r 2 dQ 2 ] 



(6) 



where f is the Eulerian radial coordinate. It is then most 
convenient to consider a single domain which is a sphere 
of radius L centered at the origin. The expressions for 
the backreaction will then correspond to this single 
domain. 

When nonlinear structures such as galaxy clusters and 
voids form during late stages of evolution of the Uni- 
verse, is the spacetime metric still perturbed FLRW? 



The answer is yes, provided matter peculiar velocities 
remain small. We showed this by considering a simple 
but generic model of spherically symmetric pressureless 
dust collapse [ll| (see also [Hj]). We now apply the back- 
reaction results © to our collapse model and show that 
the corrections are extremely small. This is the first ex- 
ample of a fully covariant calculation of the backreaction 
valid even in the nonlinear regime of structure formation. 

The spacetime of the spherically collapsing dust is 
described by the Lemaitre-Tolman-Bondi (LTB) metric 
given by 



ds 2 



-dt 2 



R' 2 dr 2 
1 — k(r)r 2 



R 2 dQ 2 . 
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Here t is the proper time measured by observers with 
i fixed coordinate r, which comoves with the dust. R(t,r) 

is the area radius of the dust shell labelled by r, and 
(3) satisfies the equation R 2 = 2GM/R - kr 2 . Here M(r) 
is the mass inside each comoving shell and a dot denotes 
a derivative with respect to t. The energy density of 
dust measured by a comoving observer satisfies p(t, r) — 
M' /4ttR 2 R' where the prime denotes a derivative with 
respect to r. Initial conditions are completely specified by 
choosing an initial density p(U,r), velocity R(U,r) and 
an initial scaling function R{U,r). In [ll| we made the 
following choices: To set the initial conditions as being a 
perturbation around the Einstein-deSitter (EdS) solution 
characterised by the scale factor a(t) = (t/to) 2 ^ 3 with 
t = 2/ (3Hq) where H is the standard Hubble constant, 
we chose the functions 



R(ti,r)=a,ir ; Rfa, r) = Oji^r , (8) 

( (1 + 5.,), r < r* 
p(ti,r)=pi< (1 — S V ), r*<r<r v (9) 
[ 1, r > r v , 

where pi = p(U), a, = (*i/*o) 2/3 , Hi = 2/(3^), and 5*, 
8 V , r* and r v are constants whose values were chosen 
so that the system being described is initially a small 
overdensity of extent < r < r*, surrounded by a small 
underdensity out to radius r v . In particular, the following 
values were chosen for the various parameters (cf. Table 

i of [HI) 



m = 0.001 
5* = 2.21 x 10~ 3 



H = 72kms~ 1 Mpc" 1 , 



S v = 5x 10 -3 . 



= 16.7 Mpc ; r v = 23.5 Mpc. 



(10) 



The exact solution for R(t, r) can be written in paramet- 
ric form in terms of trigonometric or hyperbolic functions 
depending on the sign of k(r). One then numerically ob- 
tains the function R(t, r) in the region of interest in the 
(t, r) plane. 

It was further shown in [ll[ that for the chosen model, 
for all times < t < to, the LTB metric (J7J can be 
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brought to the form 

ds 2 = -{l + 2<p)dT 2 +a 2 (T){l-2iP){dr 2 +P 2 dn 2 ), (11) 

where tp and ip satisfy \cp\, \ip\ <C 1, and a(r) = 
(3H t/2) 2 / 3 . This is achieved by the coordinate trans- 
formation (t, r) — > (r, f) given by 

r = i + £°(i,r) ; r= (fl(t, r)/a(t)) (1 + £(t, r)) . 

Here £° and £ are assumed to satisfy \£,\ "C 1 and 

are determined by integrating the equations 

C = (1/2) (k(r)r 2 + {at) 2 ) (R'/R) ; £°' = atJi? 

where € = (df/dt) m dt(R/a) is the comoving peculiar 
velocity, also assumed to remain small (13(. This is a 
self-consistent calculation, in that one assumes such func- 
tions to exist to set up the equations, and then solves for 
them showing that they satisfy the required properties. 
In particular, one finds that the peculiar velocity v does 
remain small throughout the evolution. Of course, mod- 
els with large peculiar velocities can be studied (e.g. the 
~ lGpc void studied in and such a case would indi- 
cate a breakdown of the weak field approximation. What 
is important however, is whether such inhomogencities 
are generic. The parameter choices in our model reflect 
the nature of typical inhomogeneities at the Last Scat- 
tering epoch and do not lead to large underdense regions 
and/or large peculiar velocities. With these choices, we 
find that the metric potentials <p and ip can be obtained 
self-consistently as tp — — £° + (av) 2 /2 and ip = £,°H + £, 
where, at the leading order we have H = 2/3f ~ 2/3r. 
It can be analytically shown by working in the Newto- 



nian gauge |l2j that the metric potentials tp, ip are in fact 



equal at leading order, and this can also be checked ex- 
plicitly. We emphasize that the metric potentials remain 
small in magnitude compared to unity, even though the 
density contra st ( pi p — 1) becomes completely nonlinear 
at late times [11|. While a more careful second order 
calculation can in principle be performed, it is not diffi- 
cult to show that the higher order terms will contribute 
insignificantly to the backreaction, so long as peculiar ve- 
locities are small. [This stems from the basic structure 
of the nonlinear metric potentials which depends on an 
expansion in (HR) < 1.] Hereon when considering time 
derivatives of small quantities we shall not distinguish 
between t and r. Terms quadratic in v are retained since 
on dimensional grounds one expects av ~ HR whereas 
tp, ip ~ (HR) 2 , which is confirmed by our numerics (see 
also 0). 

We now compute the backreaction given in ([3]) . These 
expressions were derived in (loj under the requirement 
that the averaging operation be free of gauge related am- 
biguities, in linear perturbation theory. However, the ac- 
tual conditions used to derive Q only depended on con- 
sidering leading order effects in the metric perturbations. 



A key step was the transformation between the metric 
(fTTI) (in Cartesian spatial coordinates) and the volume 
preserving form ((5|) , which was achieved by the transfor- 
mation t t, x A x A + d A (3 (A = 1, 2, 3), where j3 is 
the function appearing in J3J. In the present context, the 
same transformation remains valid at the leading order 
and hence the backreaction in ([3]) is physically relevant 
here as well. We emphasize that this truncated averag- 
ing operation remains valid even at late times since the 
weak field approximation for gravity works well during 
nonlinear structure formation. 

Since our numerical results involve (t, r) where r co- 
moves with the matter, we must reexpress the averaging 
operation ([4| in terms of these variables. It is easy to 
show that, at the leading order, the average of a scalar 
s(t,r) defined in (j4|) can be written as 



Ht)Lf J 



r L (t) 



sR 2 R'dr . 



(12) 



where r L (t) solves R{t,r L (t)) = a(t)L. Recall that flU} 
gives the average of s over a single domain centered at 
the origin, which is what we restrict ourselves to in this 
paper. We are constrained to consider values r < r v due 
to unphysical shell crossing singularities in the region be- 
yond (see [ill), so the largest value of L we can choose 
is L = r v , which then ensures rr,(t) < r v since r/,(t) is a 
decreasing function for this choice. This gives us an Eu- 
lerian averaging scale of L = 23.5 Mpc which is smaller 
than the more realistic expected value of ~ 100/i _1 Mpc. 
Our model potentials ip, ip and their derivatives hence 
do not strictly average to zero as needed [h], [ll[ . One 
can check that actual average values ((dAp)) 2 are small 
(< 10%) for all times compared to terms like ((dAp) 2 ) 
which are needed in the backreaction calculations, al- 
though it turns out that the time derivatives satisfy 
(ip) 2 ~ ((tp) 2 ). However, since the averaging scale chosen 
here is large enough to encompass all the inhomogene- 
ity of this system, we expect that our estimates for the 
backreaction are fairly representative. 

Consider now the function [3 which satisfies the Poisson 
equation on a flat 3-space background, such that j3 = 
if ip = = tp. We can directly write the solution for (3 in 
terms of the coordinate f as, 



I3(r,r) = -(1/4tt) / d 3 yq(T,y)/\P-y\ 



qy 2 dy 



qydy, (13) 



where q = — 2p (we have set tp = ip at leading order) and 
the integration is over Eulerian spatial coordinates. The 
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FIG. 1: The evolution of -5 (1) /QH 2 . Also shown is a 
hypothetical curvature-like correction, evolving like ~ a 



FIG. 2: The normalised evolution of the backreaction 
functions other than To enhance contrast, a strongly 

decaying early time mode for /H 2 has not been shown. 



following relations are useful in the calculations, 



in terms of t and valid at leading order, reduce to 



df/3 



i r 



qy dy, 



fig/9 = g--%3, 
r 



(14a) 
(14b) 
(14c) 



We will need o>/9, df/9 and df/3 as functions of (t,r), 
which is done by replacing r and f at leading order by t 
and R/a respectively. This gives us 



(d f l3)(t,r) = -=j / qR 2 R'dr, (15a) 
aR 2 J 

(d$f3)(t,r)=q-^ J\R 2 R'dr, (15b) 

(dim, r ) = ^[ q R2R ' dr - 2 i q+ iv ql ' (15c) 



where in the last equation we have used dfq = (a/ R')q' 
at leading order. 

Also, noting that the time derivatives in Q are taken 
keeping f fixed, we have 



(fc/3)(t,r) 



aR 2 



qR 2 R'dr, 



(16a) 



(d f P)(t,r) = — 3 / qR 2 R'dr, (16b) 
(d$p)(t,r) = q-^ J qR 2 R'dr, (16c) 



which follow from (|14[) . The expressions ([3]), rewritten 



2((^) 2 )-((a?/3) 2 )-2((l/f 2 )(^/3) 2 ) 



6<(5^) 2 ) - (($/9) a ) - 6<(fi ? /3 - fdUf/f^ 



77(1) + t?(2) = _ 2 ff 



((5?/9)(t??/9)) + 2<(l/f 2 )(fi r -/3)(d r -/9)> 



£( 2 ) = ((5 f /3 + Hd?l3)(a 2 Hd?l3 - Bpip)) , 



(17) 



where the angular brackets are now defined by (TT2")) and 
the various integrands can be read off using (fl~5"j). P^|) 
and the results df<p « (a/R')ip' and f « (R/a). 

Figs. [T] and show results of numerical calculations 
performed with Mathematica. Fig. [T] shows the evolu- 
tion of the dominant correction —S^/QH 2 , as a func- 
tion of the scale factor. The dashed line shows a hypo- 
thetical curvature like correction. The actual backreac- 
tion evolves differently due to significant evolution of ip. 
Note that the largest value of \S^/H 2 \ computed here 
is ~ 10~ 6 , whereas estimates using linear theory suggest 
this value should be ~ 10 -4 [ioj . This discrepancy high- 
lights an issue noted in [Io| , namely that nonlinear inho- 
mogeneities on small scales do not contribute significantly 
to the backreaction. Our model has no large scale inho- 
mogeneities and underestimates the backreaction. Reas- 
suringly, accounting for the deficit only requires a calcu- 
lation in linear theory, such as the one in (loj . Fig. [2] 
shows the evolution of the remaining (normalised) inte- 
grals. An initial rapid decay of V^/H 2 starting from 
~ 10~ 8 has not been shown, in order to enhance the con- 
trast in the late time behaviour of the three functions. 
The other functions remain subdominant compared to 
at the early times not shown. 

Our covariant and self-consistent calculation of the 
backreaction in this spherical collapse model establishes 
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that inhomogeneities have an insignificant impact on the 
average cosmological dynamics. In particular, the ob- 
served cosmic acceleration cannot be explained by the 
averaging of inhomogeneities. Our nonlinear dust model 
can be regarded as representing a realistic situation, be- 
cause it has a overdensity-void structure, and departure 
from sphericity, tidal interactions, and second order cor- 
rections are not expected to introduce any significant 
change in the results. What appears true in general is 
that as long as peculiar velocities remain small, as seems 
to be the case in the real Universe, a description as a 
perturbed FLRW model is valid, and this keeps the back- 
reaction small. 
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